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The turbulent transport of minority species/impurities is investigated in 2D 
drift-wave turbulence as well as in 3D toroidal drift- Alfven edge turbulence. The 
full effects of perpendicular and - in 3D - parallel advection are kept for the 
impurity species. Anomalous pinch effects are recovered and explained in terms 
of Turbulent EquiPartition (TEP) 



1 Anomalous Pinch in 2D Drift- Wave Turbulence 

The Hasegawa-Wakatani model Q for 2D resistive drift-wave turbulence reads 

d t (n-x) = c(cp-n)+^,Vj_n, d t (0 = C(q>-n) + fi (a V 2 1 _(o, (2) 

with co = V^cp and d t = d/dt + z x Vj_cp ■ Vj_. Here, n and cp denote fluctuations 
in density and electrostatic potential, co is the vorticity, V x z x Vj_(p. The pa- 
rameters in the HW system are the parallel coupling C, and diffusivities, /u n ^a>- 
2D impurity transport in magnetized plasma is modeled by the transport of a pas- 
sive scalar field: 

d t Q - • (e^V ± <pJ = ^eVie, (3) 

where is the density of impurities, /jq the collisional diffusivity, and C, = ^-^jf- 
the influence of inertia, which enters via the polarization drift. The latter makes 
the flow compressible, consequently for ideal (massless) impurities, C, = and 
advection is due to the incompressible electric drift only. In all cases the impurity 
density is assumed to be so low compared to the bulk plasma density that there is 
no back-reaction on the bulk plasma dynamics. 
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Figure 1 : (a) Vorticity and (b) density of inertial impurities in the saturated state 
with C = 1 and £ = 0.01, L = 40. Other parameters: /u n = ^co = A*e = 0.02. 



1.1 Vorticity - Impurity correlation 

The equation for the impurities can be rewritten in the form: 

J r (ln0-^co) = CV ± ln0-J f Vj_cp + ^v2e. 

9 

If the diffusivity /uq is of order £ <C 1 and fluctuations 9i of the impurity density 
measured relative to a constant impurity background 9o do not exceed a corre- 
sponding level, the quantity In 9 — C,oo is approximately a Lagrangian invariant. 
Turbulent mixing will homogenize Lagrangian invariants in TEP states GlEl, 
leading to 

In 9 — ^(0 ~ const, 

which constitutes a prediction about the effect of compressibility on the initially 
homogeneous impurity density field. The conservation of impurity density yields 

1 r 

— «i + Go>, 

which conforms with the assumed ordering. We thus predict a linear relation be- 
tween impurity density 9 and vorticity (0, the proportionality constant being the 
mass-charge ratio This is related, but not the same as, to the aggregation of 
dense particles in vortices in fluids due to the Coriolis force |H). The predic- 
tion is verified by numerical simulations of inertial impurities in saturated HW- 
turbulence for C = 1. The simulations are performed on a [—20, 20] 2 domain, 
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Figure 2: Scatter plot of impurity 
density and the vorticity field at t = 
100 for different values of the mass- 
charge ratio C, in the saturated state in 
HW with c = 1: £ = 0.05 (red), £ = 
0.01 (green), and £ = 0.002 (blue). 
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Figure 3: Evolution of the radial 
drift velocity of inertial impurities in 
the saturated state in HW with C = 
1 . The impurities are uniformly dis- 
tributed at t = 0. 



using 512 2 gridpoints, and impurity diffusivity 0.02. The impurity density field is 
initially set to unity. The impurity density field for C, = 0.01 is presented together 
with vorticity in Figure U Figure |2] shows a scatter plot of the point values of 
impurity density and vorticity at time 150 for three different values of The pro- 
portionality factor = 1 + Ku> is determined to be slightly below one: K ~ 0.82^. 



1.2 Anomalous pinch 

The role of inertia for a radially inward pinch is investigated by considering the 
collective drift of impurities. Ideal impurities do on average not experience a drift, 
but this is not the case for inertial impurities, since compressibility effects arrange 
for a correlation between 0i and (0. Note that only the deviations from the above 
discussed linear relationship = 1 +K(D result in a net flow, as / K(tiv r dx = for 
periodic boundary conditions. 

The evolution of the radial drift velocity, measured as the net radial impurity 
transport, is presented in Figure[3] The radial drift velocity has a definite sign that 
depends on the sign of There is a continuous flow of impurities in a definite di- 
rection (inward for positively charged impurities). This resembles the anomalous 
pinch observed in magnetic confinement experiments [5 1. Average radial drift ve- 
locities computed using the values of the drift from t — 25 to t = 150 are presented 
in Table [T] The scaling of the average radial drift with C, is seen to be remarkably 
linear. 
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Table 1 : Radial drift velocity of impurities for different values of the mass-charge 

ratio £ in the saturated state in HW with C = 1. Calculated as the average value 

between t = 25 and t = 150. Parameters: /u n — fi a = /jq = 0.02. 
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2 Drift- Alfven Turbulence 

We now consider drift- Alfven turbulence in flux tube geometry J51|71[H1. The 
following equations for the fluctuations in density n, potential (b with associated 
vorticity co = V^_(b, current J and parallel ion velocity u arise in the usual drift- 
scaling: 



+ {§, co} = (n) + V„7 -f^coV^o), 



-^- + {§,n EQ + n} = % (n-cb)+V n (/-M)+// n V^n, 
^ (pM, , + Jif) + $ {(b, J} = V, , (n EQ + n - (b) - CJ, 
e(^ + {<M}^) =-V„(n £G +n). 
The evolution of the impurity density is given by 

d t Q = (£/e)v ± . (e^v ± (b) - e^: ((b) - v., (e«) - ^vie 



(4a) 
(4b) 
(4c) 
(4d) 

(5) 



Standard parameters for simulation runs were Ji = 5, q = 3, magnetic shear s~= 1, 
and cog = 0.05, with /u w = [i n = 0.025, corresponding to typical edge parameters 
of large fusion devices. Simulations were performed on a grid with 128 x 512 x 32 
points and dimensions 64 x 256 x 2% in x,y,s corresponding to a typical approx- 
imate dimensional size of 2.5 cm x 10 cmx 30 m [6|. Here we present results 
from a low j3 = 0. 1 run with C = 11.5. In FigureH|the dynamical evolution of the 
impurity density is exemplified in a plot showing the poloidal projection of the 
impurity density. 
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Figure 4: Impurity distribution projected onto a poloidal cross-section (radial di- 
mension not to scale). Initial distribution (left) and after 25 time units (right). 

The flux r of the impurity ion species can in lowest order be expressed by the 
standard parameters used in modeling and in evaluation of transport experiments: 
a diffusion coefficient D and a velocity V, which is associated to a pinch effect, 

r y (s) = -D( S )d x {Q) y + V( S ){Q) y . (6) 

From scatter plots of T(r)/ (n) y versus d x ln(n) y , values for D(s) and V(s) are 
obtained. The poloidal (coordinate s) dependence of D and V is rather strong and 
shown, with numerical uncertainties, in Fig. The effective advective velocity 
V(s) changes sign and is at the high field side directed outwards. This pinching 
velocity is due to curvature and can be consistently explained in the framework 
of Turbulent EquiPartition (TEP) @ |U as follows: In the absence of parallel 
advection, finite mass effects and diffusion, Eq. © has the following approximate 
Lagrangian invariant 

L(s) =\nQ + a> B xcos(s) — (£>Bys'm(s) . (7) 

TEP assumes the spatial homogenization of L by the turbulence. As parallel 
transport is weak, each drift plane s = const, homogenizes independently. This 
leads to profiles (L(s)) y = const. (s). At the outboard midplane (s = 0) the im- 
purites are effectively advected radially inward leading to an impurity profile 
((lnG)^ oc const. — (Dgx), while at the high field side they are effectively advected 
outward ((InG) const. + QOgx). One should note that this effective inward or 
outward advection is not found as an average E x B velocity, but is mitigated by 
the effect of spatial homogenization of L under the action of the turbulence. The 
strength of the "pinch" effect is consequently proportional to the mixing properties 
of the turbulence and scales with the measured effective turbulent diffusivity. We 
arrive at the following expression for the connection between pinch and diffusion: 

V(s) = -au) B cos(s)D(s) . (8) 
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Figure 5: Impurity diffusion D (a) and pinch velocity V (b) over poloidal position 
(s) with error-bars. The pinch velocity is compared to QOb * cos(s) * D(s) (dashed 
line). 



Considering a stationary case with zero flux and Eq. © we obtain a = 1. The bal- 
looning in the turbulence level causes the inward flow on the outboard midplane 
to be stronger than the effective outflow on the high-field side. Therefore, aver- 
aged over a flux surface and assuming a poloidally constant impurity density, a net 
impurity inflow results. This net pinch is proportional to the diffusion coefficient 
D in agreement with experimental observations ifTOl . 
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